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Effect of nearby supernova remnants on local Cosmic-Rays 
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ABSTRACT 

We study in detail the effect of different particle release times from sources on the cosmic -ray 
(CR) spectrum below 10 15 eV in the Galaxy. We discuss different possible forms of particle 
injection such as burst-like injection, continuous injection for a finite time, injection from a 
stationary source and energy dependent injection. When applied to the nearby known super- 
nova remnants, we find that the observed CR anisotropy data favour the burst-like particle 
injection model for the CR diffusion coefficient D(E) oc E a with a = 0.3 — 0.6 in the local 
region. In the study we have also found that the contribution of the sources Gl 14.3+0.3 and 
Monogem dominate if the observed anisotropy is a result of the effect of the nearby sources. 
Further study shows that we should not neglect the contribution of the undetected old sources 
to the local CR anisotropy. 
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1 INTRODUCTION 



Considering the similarity of the power supplied by a supernova 
explosion (SNE), which is ~ 10 42 erg/ s, and the power required 
to maintain the cosmic-ray (CR) energy density in the Galaxy 
(~ 10 41 erg / s), it is widely beleived that the majority of CRs in the 
Galaxy are accelerated in supernova remnant (SNR) shock waves. 
If this is true, then it is very important to know the maximum en- 
ergy of CRs that can be accelerated in SNRs as well as the fraction 
/ of SNE energy that can be converted into CRs by diffusive shock 
acceleration. Studies based on the plane-wave approximation have 
shown that the maximum energy can be~Zx 10 14 eV, where Z 
is the charge of the particle (see e.g. Lagage & Cesarsky 1983). 
However, studies using nonlinear effects produced by accelerated 
CRs have found that the maximum energy can be as large as 
~Zx 10 15 eV r (Berezhko 1996). Regarding the total energy trans- 
fered into CRs in a SNE, there can actually be a considerable vari- 
ation in the value of / among the sources, but its value is generally 
considered to be ~ 0.1 which can easily be up to ~ 0.3 (Berezhko 
et al. 1996). 

An interesting feature of the CRs observed at the Earth is 
that they are highly isotropic. Their anisotropy amplitude is only 
~ (10~ 4 - 10" 3 ) in the energy range of (10 11 - 10 15 )eV (Guillian 
et al. 2007 and references therein) with the phase mainly found in 
the outer Galaxy particularly in the second quadrant of the Galaxy. 
The high degree of isotropy is because CRs do not reach Earth 
from the sources directly. They undergo difusion in the Galaxy as 
a result of scattering either by magnetic field irregularities or by 
self-excited Alfven and hydromagnetic waves before reaching the 
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Earth. In the present study, because CRs are composed mainly of 
protons below the knee region (~ 3 x 10 15 eV), the anisotropy 
below the knee is considered mainly because of the proton compo- 
nent. 

The anisotropy observed at the Earth may be because of the 
effect of global diffusion leakage of CRs from the Galaxy, the effect 
of random nature of CR sources both in space and time, or the effect 
of the local sources (see e.g. Ptuskin et. al. 2005). Studies using 
local SNRs lying within lkpc from us with ages less than 10 yr 
have been carried out, assuming these to be instantaneous point 
sources (Dorman et al. 1984). It was concluded that the fact that the 
CR anisotropy is almost energy independent in the energy range 
of (10 12 — W 14 )eV may be because we are observing the Vela 
source's maximum anisotropy amplitude. 

In a recent paper, Erlykin & Wolfendale (2006) have shown 
that the observed characteristics of CRs, and in particular their 
anisotropy, are mainly the product of local sources located within 
a distance of (1 — 2)kpc from the Sun. Their study considers the 
actual positions of the observed SNRs within 3fcpc from the Sun 
as well as the position of the solar system at the inner edge of the 
Orion arm. However, they could not incorporate the actual ages of 
the SNRs in their study because their model considers the CRs to 
be confined in the SNRs for around 8 x W 4 yr, but the majority 
of observed SNRs have ages less than ~ 4 x 10 4 yr. Among the 
ten samples that they simulated, only one sample has been found to 
show the phase of the anisotropy which agrees with the observed 
data. 

Considering local SNRs located within a distance of lkpc and 
age less than 5 x W 4 yr, a study had also been carried out using the 
data on the high energy CR spectrum and anisotropy to constrain 
the value of CR diffusion coefficient in the Galaxy and its depen- 
dence on the energy (Ptuskin et al. 2005), assuming the particles 
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are released from the sources at the time of the explosion itself. 
Ptuskin et al. have shown that the diffusion model with reaccelera- 
tion is barely compatible with the observed CR anisotropy, with a 
discrepancy that is roughly within a factor of 3. The plain diffusion 
model predicts a value of anisotropy that is too large. 

From these earlier works, it can be seen that the cause of 
the observed amplitude and the phase of the CR anisotropy is not 
clearly understood, although local phenomena are more likely to 
play a dominant role. In the present study we also try to under- 
stand the problem, but with more emphasis given to the particle 
release time. We consider different possible forms of particle injec- 
tion from the sources, such as like burst-like injection, continuous 
injection for a finite time, injection from a stationary source and en- 
ergy dependent injection. We then study the observed local SNRs 
and examine the effect of the different particle release times on the 
observed CR anisotropy. In short, the plan of the paper can be out- 
lined as follows. In Section 2 we calculate the CR spectrum gener- 
ated by a single point-like source. In Section 3 we discuss the effect 
of different particle release times from the sources on the CR flux 
and anisotropy. In Section 4 we give an application to the nearby 
SNRs located within 1.5kpc from the Earth. Finally, in Section 5 
we present the results and discussions about the whole study. 



2 CRS FROM A SINGLE POINT LIKE SOURCE 
2.1 CR spectrum 

In the diffusion model, neglecting convection, energy losses and 
particle losses as a result of nuclear interactions, the propagation 
of CR protons in the Galaxy is given by (see e.g. Gaisser 1990 and 
references therein) 

ON 
~dt 

where N(r, E, t) is the differential number density, E is the pro- 
ton kinetic energy, D(E) oc E a with a — constant (positive) 
is the diffusion coefficient which is assumed to be spatially uni- 
form in the Galaxy and Q(r,E,i) is the proton production rate 
[i.e. Q(r, E, t)d 3 rdEdt is the number of protons produced by the 
source in a volume element d 3 r in the energy range (E, E + dE) 
in time dt]. 

The Green function G(r, r', t, t') of Eq. (1) [i.e. the solution 
for a (5-function source term Q(r,t) = 8(r—r')8(t—t')] is obtained 
as 



V • (DVJV) +Q = 



(1) 



G(r,r',M') = 



871- 3 / 2 [D(t - t')] 



:-s/2 



exp 



4D(t - t') 



Then, the general solution of Eq. (1) is given by 

/OO ft 
dr dt'G(r,r',t,t')Q(r',E',t') 
-co J — CO 



(2) 



(3) 



Now, to calculate the number density N(r, E, i) at a distance r 
from a point source, let us consider a small source volume ? where 
protons are generated. Before the protons are released into the inter- 
stellar medium (ISM), they are assumed to be distributed uniformly 
in the source volume with density qo(E). In other words, we can 
write 



Q(r,E,t) = q (E)q(t), 
= 0, 

and for a point source 



inside the source, 
outside the source, 



(4) 



lim ? _ Qo{E)<; = q(E) 

With this assumption the proton density at a distance r from a 
point source can be calculated by introducing Eqs. (2)&(4) in Eq. 
(3) as 



N(r,E,t) 



8tt 3 / 
xexp 



- 3 / 2 y_< 



dt' 



[D(t - t')f /2 



4D(t - t') 



q{E)q{t') 



(5) 



and the proton flux is given by 1(E) « (c/4ir)N(E), where c is 
the velocity of light. The proton source spectrum q(E) is taken as 



q{E) = k(E 2 + 2Em p y {r+1)/2 (E + rn p ) 



(6) 



where ft is the normalization constant and m p is the proton mass 
energy. The source spectral index V is chosen such that F + a — 
2.73, the observed proton spectral index (Haino et al. 2004). In Eq. 
(5), the function q(t') carries the information about how and when 
the CR particles are emitted by the source. The actual form of q(t') 
is not known exactly. It can be a simple energy independent func- 
tion or a complicated energy dependent function also (i.e. particles 
with different energies emitted at different times). In fact, detailed 
studies have shown that the most energetic paricles start to escape 
from the source region at the very beginning of the Sedov phase 
itself, although the major fraction of the CRs remain confined for a 
period of about 10 5 yr (Berezhko et al. 1996). In the present study, 
we will consider different possible forms of q(t')'s corresponding 
to different types of particle injection, that is, burst-like injection, 
continuous injection for a finite time, injection from a stationary 
source and also the energy dependent injection mentioned above. 



2.2 CR anisotropy 

The amplitude of CR anisotropy due to a single source in the diffu- 
sion approximation is given by (Mao & Shen 1972) 



5i = 



3D |VjV» 
c Ni 



(7) 



where Nj is given by Eq. (5) for a point source i with distance r ; 
and age U. The total anisotropy parameter at the Earth, due to a 
number of nearby discrete sources in the presence of an isotropic 
CR background, is given by 



^ Ii8jfj.hr, 



It 



(8) 



where the summation is over the nearby discrete sources. 
fi denotes the direction of the source i giving a flux E 
and fm denotes the direction of maximum intensity. It = 
1.37(E/GeV)- 2 - r:i cm- 2 s- 1 sr- 1 GeV- 1 represents the total 
observed flux of CR protons above ~ WGeV (Haino et al. 2004). 
The direction (phase) of the anisotropy is given by the direction 
of maximum intensity. Therefore, the anisotropy 6 as well as the 
phase at an energy E depends on the age and distance of the nearby 
sources. These can be determined by different sources at different 
energy intervals. In fact, the energy dependence of 8 depends on 
the dependence of Ii, I T and 5i on energy and are discussed in de- 
tail in the next section for different particle injection forms. In case 
a single source is dominant over the entire energy range, the total 
anisotropy S is given by 



r m r 

= —Or, 

It 



(9) 
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where the subscript m denotes the source giving the maximum flux 
at the Earth. 



3 EFFECT OF PARTICLE RELEASE TIME FROM THE 
SOURCES 

3.1 Burst- like injection 

When the injection time of CRs from the sources is significantly 
short, the injection process can be represented by a 5-function as 



q(t') = S(t' - to) 



(10) 



Eq. (10) represents the burst-like injection of all the particles at time 
to after the explosion. Using this, the proton density is obtained 
from Eq. (5) as 



N(r,E,t) 



q b {E) 



8tt3/ 2 [D(t - t )] 3/2 



exp 



4D(i- to) 



(11) 



For the present case, the source spectrum qb(E) is given by Eq. 
(6) with the constant k obtained by normalizing the proton source 
spectrum such that the total amount of kinetic energy contained in 
the injected protons is / ~ (0.1 — 0.3) times the total SNE energy. 
With this, the single source anisotropy (Eq. 7) becomes 



Si 



2c (U - to) 



(12) 



which is independent of energy. For the high energy particles for 
which the diffusive time td ~ r 2 / D taken to travel a distance r 
in the Galaxy is less than the time interval after their release from 
the source (t — to), i.e. td < [t — to), the exponential term in Eq. 
(11) tends to 1 and hence, the proton spectrum follows a simple 
power law as N(E) oc £:~( r + 3a / 2 ). Similar result can also be 
found in Aharonian & Atoyan 1996 where they studied about the 
pion decay gamma radiation in the vicinity of the CR sources. Also, 
the total proton spectrum observed at the Earth It can be written in 
the form I T (E) oc E~ < - r+a) where V + a = 2.73 as mentioned in 
section 2.1. Therefore, the total anisotropy 5 (Eq. 9), when a single 
source dominates, is found to follow an energy dependence of the 
form 8 oc E~ a ^ 2 . S also depends strongly on the time to after the 
SNE at which particles are injected into the ISM. Applying Eqs. 
(1 1)&(12) into Eq. (9), the dependence of high energy particles on 
injection time is given by S oc (t m — to)~ 5 ^ 2 where t m denotes the 
age of the source giving the maximum intensity at an energy E at 
the Earth. This shows that as the value of to increases, the total CR 
anisotropy increases. 



3.2 Continuous injection for a finite time 

For the case of continuous injection of CRs from a source for a 
finite time interval from ti to ti where (0 ^ ti < ti), the source 
is switched on only between t\ and ti. So, the temporal particle 
injection function can be written as 



q(t') = 1, for ti ^ t! ^ ti, 
= 0, otherwise, 

The source spectrum (Eq. 6) for this case is given by 

q b (E) 



q c {E) = 



(t 2 -ti) 



(13) 



(14) 



where qt(E) has been obtained in section 3.1. The proton density 
can then be easily obtained from Eq. (5) as 

N(r,E,t) = ^[erf(V^)-erf(V^)] (15) 



where, 



and X2 — -777 r , (16) 



4D(t-ti) " ' AD(t-t f ) 

er f(y) — (2/^/7?) J* exp(—u 2 )du is the error function and tf = 
min[t,t2]. The corresponding single source anisotropy amplitude 
(Eq. 7) is given by 



3D 
Oi = — 



fx\exp(—x\) — ^/xiexp(-xi) 
erf(V^) ~ er f(V^i) 



(17) 

It is possible to discuss the particle spectrum given by Eq. (15) for 
two seperate cases: ti < t < ti and ti < ti < t. The two cases 
are expected to give completely different spectra, particularly for 
the higher energy particles as discussed below. 

Ca.se 1: (ti < t ^ ti) : — For this case, we observe the 
source while it is still active and liberating CRs continuously into 
the ISM. Here, tf = t which implies X2 — oo (using Eq. 16). 
By the property of error function, erf(^/xi) — 1 for X2 = oo 
(Abramowitz & Stegun 1964). So, Eq. (15) for this case can be 
written as 

Qc(E) 



N(r,E,t) 



[1 - er/(V5D] 



(18) 



A-KrD 

For the high energy particles with a diffusion time t d « (t — ti), 
x\ — » and because erf(^/xi) —* for x\ — > 0, the 
spectrum given by Eq. (18) follows a power-law of the form 
N(E) oc _E~( r+a ) and the single source anisotropy (Eq. 17) 
has an energy dependence given by Si oc E a . Thus, the total 
anisotropy (Eq. 9) varies with energy as 8 oc E a in the high energy 
regime in this case. 

Case 2: (ti < ti ^ t) : — For this case, tf = ti and 
hence, X2 = r 2 /4D(t — ti). Here, the proton spectrum is given by 
Eq.(15), but with tf replaced by ti. For very high energy particles 
for which x\ — * 0, we can safely write xi « 1 as t\ < ti. How- 
ever, a property of the error function gives erf(^/xi) w 2^xi/n 
for y/xi << 1. So, for high energy particles the spectrum can be 
written as 



N(r,E,t) 



Qc(E) 



(19) 



(20) 



4(7rD) 3 / 2 Vt T_ ^ 
and the single source anisotropy can be written as 

S ~ 3 r% 
1 ~ 4c [U - t 2 ) 

The energy dependence of Eqs. (19)&(20) is the same as in the case 
of burst-like injection of particles discussed in section 3.1. This 
implies that the total anisotropy parameter (Eq. 9) varies as 5 oc 

E~ a/2 . 

From the study, we can now see that for the continuous in- 
jection of particles it is important to know when the observation 
is being made for a particular source. If we observe the source af- 
ter it is switched off, we see a spectrum much steeper than would 
have been observed during its active period. For the present study, 
we take ti = 10 5 yr because it is generally accepted that CRs re- 
main confined in the SNRs upto around 10 s yr (see section 3.4). 
So, for the nearby SNRs which are listed in Table 1, it is expected 
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that for the same source spectrum the Geminga supernova will give 
a steeper spectrum at the Earth compared to all the other remain- 
ing SNRs, because it is the only source (with an estimated age of 
t ~ 3.4 x 10 5 yr) that is being observed after it is switched off. 
Another important point is that for a fixed value of t2, the total 
anisotropy parameter 5 (Eq. 9) varies with ti. In both cases dis- 
cussed above, this is given by S oc l/(fa — it) for high energy 
particles. Thus, the dependence of S on injection time in the con- 
tinuous particle injection model is weaker than that in the burst-like 
injection model. 



3.3 Injection from a stationary source 

The particle spectrum for the case of injection of CR particles from 
a stationary source can be easily obtained using Eq. (11). For the 
burst-like injection of particles at time to = 0, Eq. (11) gives 



N(r,E,t) = 



q b {E) 



8n 3 / 2 (Dt) 3 / 2 



exp 



iDt 



(21) 



For a stationary source injecting a q s (E) number of particles of en- 
ergy E per unit time continuously, the particle density N(r, E) is 
obtained by replacing qt(E) in Eq. (21) with q s (E)dt and integrat- 
ing over dt from to oo as 



N(r,E) = 



qs(E) 
4%rD 



(22) 



This shows that the particle spectrum follows a simple power-law 
of the form N(E) oc _E~( r+a ) as in the continuous injection case 
(Case 1 of section 3.2). The single source anisotropy amplitude in 
this case is 



St = 



3D 



(23) 



and the total anisotropy (Eq. 9) follows S oc E a . 



3.4 Energy dependent particle injection 

Studies based on diffusive shock acceleration in SNRs have shown 
that the highest energy particles already start leaving the source re- 
gion at the beginning of the Sedov phase (Berezhko et al. 1996), but 
the major fraction of accelerated CRs remain confined for almost 
around 10 5 yr for an ISM hydrogen atom density of nu = 1cm -3 . 
Thus, an energy dependent proton confinement time in the source 
region (proton release time from the source) can be written as 
(Berezhko & Volk 2000) 



T con {E) = mm[10 5 , 10 3 (E / E max y 5 ]yr 



(24) 



where E max , the maximum energy of accelerated particles, is taken 
as 10 15 eV (Berezhko 1996) and 10 3 yr marks the start of the Sedov 
phase. For such a type of particle injection, the function q(t') in Eq. 
(5) can be written as 



q(t') = S(t' — Tcon) 



(25) 



The particle spectrum and the anisotropy parameter for this case 
are the same as in the burst-like injection case with to replaced by 
T co „(E). From Eq. (24), it can be seen that all the particles with 
energies less than ~4x 10 14 eV remain confined until 10 yr. So, 
except for particles with £ ) 4x 10 14 eV both the spectrum and 
the anisotropy are the same as those of the burst-like injection of 
particles at time to = W 5 yr. 



Table 1. Parameters of SNRs with known ages located within a distance of 
1.5kpc from the Earth (References: 1 Braun et al. 1989; 2 Strom 1994 and 
for other sources, see the references given in Thoudam 2006). 



SNR 


Kdea ) 


d(hpc) 


t(iiv\ 


G65.3+5.7 


65.3 


1.0 


14000 


G73.9+0.9 


73.9 


1.3 


10000 


Cygnus Loop 


74.0 


0.4 


14000 


HB21 


89.0 


0.8 


19000 


Gl 14.3+0.3 


114.3 


0.7 


41000 


CTA1 


119.5 


1.4 


24500 


HB9 


160.9 


1.0 


7700 


S147^ 


180.0 


0.8 


4600 


Vela" 


263.9 


0.3 


11000 


G299.2-2.9 


299.2 


0.5 


5000 


SN185^ 


315.4 


0.95 


1800 


Monogem 


201.1 


0.3 


86000 


Geminga 


195.1 


0.4 


340000 



4 APPLICATION TO NEARBY SNRS 

The diffusion coefficient D(E) oc E a , which governs the diffusion 
of CRs in the Galaxy, is generally obtained using the observed sec- 
ondary to primary boron-to-carbon (B/C) ratio. However, the values 
of D(E) obtained from the same experimental data are different for 
different CR propagation models in the Galaxy. The values vary 
from D(E) = 2 x W 28 (E /5GeV) 0S cm 2 s' 1 for E > 5GeV, 
where E is in GeV (Engelmann et al. 1990), obtained using a leaky 
box propagation model, to D(E) = 5.9 x 10 28 /3E C 13 cm 2 s" 1 , 
where f3 = v/c and v is the velocity of the particle (Jones et 
al. 2001), given by a diffusion model with particle reaccelaration. 
Here, we consider both cases for the present study with the source 
spectral index F chosen such that F + a — 2.73. 

Considering SNE at various distances from the Earth, a recent 
study has shown that only those sources located within a distance 
of ~ 1.5kpc can produce appreciable temporal fluctuations in the 
CR proton flux observed at the Earth (Thoudam 2006). So, for the 
present study we consider only those SNRs whose ages are known 
and which are lying within l.hkpc from us. These SNRs are listed 
in Table 1 along with their galactic longitudes (I), distances (d) 
and ages (t). The SNR RXJ1713.7-3946 is not included in the list 
because of the uncertainty in its distance. Its distance of (6± l)kpc 
is in contrast to the distance of lkpc estimated from the soft X-ray 
obsorption (Koyama et al. 1997). 

Among the various particle injection schemes considered in 
section 3, only the burst-like and the continuous injections are con- 
sidered for application to the nearby SNRs. The case of injection 
from a stationary source is not considered, mainly because every 
source has a finite lifetime and should also have a finite particle 
injection time. Moreover, the proton spectrum and the anisotropy 

5 follow the same energy dependence as in the case of continu- 
ous injection for a finite time period (see section 3.3). Similarly, 
energy dependent particle injection is the same as burst-like in- 
jection at time to — 10 5 yr except for particles with energies 
E ^ 4 x W 14 eV as discussed in section 3.4. 

Figs. 1&2 show the comparision of the CR anisotropy ampli- 
tude observed experimentally with the contribution of the SNRs 
listed in Table 1 calculated using the continuous and burst-like 
particle injection models, respectively. The experimental data are 
taken from the compilation of various results given in Erlykin & 
Wolfendale 2006. The model calculations use Eq. 8 and assume 
D{E) oc E a with a = 0.3 or 0.6 and a total SNE energy of 
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Kinetic energy E (GeV) 




1000 10000 100000 1e+06 

Kinetic energy E (GeV) 

Figure 1. CR anisotropy <5 (Eq. 8) for the continuous injection case. The 
solid lines represent the anisotropics calculated for ti = (0, 10 3 ,5 X 
10 3 , 1 X 10 4 , 2 X 10 4 , . . . . , 9 X 10 4 )yr at a fixed value of t 2 = W 5 yr. The 
fraction of the total SNE energy of 10 51 ergs going into the kinetic energy 
of protons is taken as / = 0.1. Data points are taken from the compilation 
of different experiments given in Erlykin & Wolfendale 2006. Top : For a 
diffusion coefficient D{E) oc E a with a = 0.3. Bottom : a = 0.6. 



10 51 ergs with a fraction / = 0.1 of it going into the kinetic en- 
ergy of the injected protons. From Fig. 1 (bottom), it can be seen 
that for a = 0.6 the model with continuous injection of particles 
gives an anisotropy that is too large to explain the observed value. 
However, for a — 0.3 and ti < W 4 yr, the continuous injection 
model predicts a value which agrees with the data within a factor 
of 2 for energies below 10 13 eV and a much larger value above it, as 
shown in Fig. 1 (top). However, Fig. 2 shows that the burst-like par- 
ticle injection model is able to explain the data more closely, both 
at a = 0.3 and 0.6 for some values of to. For a — 0.3, considering 
the acceptable range of / = (0.1 — 0.3), it is found that the data are 
explained at almost two different values of to: to ~ 2 x 10 4 yr and 
to ~ 5 x I0 4 yr. For to = 2 x W 4 yr, the total anisotropy 5 is deter- 
mined by the Monogem SNR for energies E < 10 13 el/andby the 
SNR Gl 14.3+0.3 for E > 10 13 eV, as shown by the dot-dashed 
line in Fig. 2 (top). This also results in a change of the direction 
(phase) of 5 at E ~ 10 13 eV\ It should be noted that the position 
of the Monogem SNR in the Galaxy with longitude I — 201.1° as 
well as that of the SNR Gl 14.3+0.3 (/ = 114.3°) are close to the 
direction of the observed phase of the anisotropy (the second quad- 
rant of the Galaxy). However, for to ~ 5 x 10 4 j/r, the data are ex- 
plained solely by the Monogem SNR, as shown by the three dotted 
lines in Fig. 2 (top). For the other case (i.e. a — 0.6), Fig. 2 (bot- 
tom) shows that the burst-like injection model with to ~ 5 x 10 4 yr 
explains the data effectively above 10 12 eV while below 10 12 eV 
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Figure 2. CR anisotropy 5 (Eq. 8) for the burst-like injection case. The lines 
represent the anisotropies calculated at to = (0, 10 3 , 5 x 10 3 , 1 X 10 4 , 2 X 
10 4 , ....,10 X W 4 )yr. Other model parameters and the data points are the 
same as in Fig. 1. Top : Diffusion coefficient D(E) oc E a with a = 0.3. 
The three dotted lines calculated at to = 4 X 10 4 yr (bottom-dotted line), 
5 X 10 4 yr (middle-dotted line) and 6 X 10 4 j/r (top-dotted line) and the 
dot-dashed line (calculated at 2 X 10 4 yr) are those which closely explain 
the observed data for the acceptable range of / = (0.1 — 0.3) (see text 
for details). Bottom : a = 0.6. The dotted line represents the anisotropy 
calculated at to = 5 X 10 4 yr which explain the data most effectively in 
this case. 



the model predicts a value that is larger than the observed value by 
a factor of about 4 (as shown by the dotted line). In this case, the 
anisotropy is also determined by the Monogem SNR, as in the case 
of a — 0.3 and to ~ 5 x lQ 4 yr discussed above. 

The figures, particularly Fig. 2, show that depending on the 
age and distance of the local sources, the total anisotropy can be 
determined by different sources at different energy ranges. These 
can be seen as breaks in few of the lines. Also note the similarity of 
the lowest lines in the bottom plots of Figs. 1 & 2. AlThough the 
two lines correspond to two completely different particle injection 
models, the similarity is because, for a continuous injection model 
with ti = 9 x lQ 4 yr (and ti = 10 5 yr) and a burst-like injection 
model with to = W 5 yr, only the Geminga supernova with an es- 
timated age of t ~ 3.4 x 10 a yr has liberated the CR particles in 
the local region (see Table 1). Moreover, as already mentioned, the 
high energy particle spectrum and the anisotropy parameter S from 
Geminga in this case of continuous injection are the same as in the 
burst-like injection case (see Case 2 of section 3.2). 

Another probable contribution to the local CR anisotropy that 
we have not considered yet is that from the undetected old SNRs. 
Studies assuming an adiabatic phase in the SNR evolution have 
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shown that the surface brightness estimated from an SNR of age 
~ 10 5 yrs lies below the detection limit of radio telescopes (Ko- 
daira 1974, Leahy & Xinji 1989). Such old SNRs can contribute 
a significant flux of CR protons at the Earth, although the elec- 
trons present inside it are not energetic enough to produce a strong 
radio flux that can be detected by present day telescopes. Fig. 3 
shows the anisotropies due to such old sources calculated by tak- 
ing their ages t = 2 x 10 5 yr (below which they are assumed to 
be detectable) and a = 0.3 or 0.6 assuming burst-like injection of 
particles at to = 10 yr. It is worth mentioning at this point that 
for such old sources the particle injection stops long before the ob- 
servation time because CRs are assumed to remain confined in the 
sources only upto around W 5 yr in the present model, and hence the 
particle injection can be well approximated by the burst-like model 
(see Case 2 of section 3.2). The anisotropy at each energy point in 
Fig. 3 is calculated using Eq. 8 and by taking the source position 
at a distance r = r max (E), where r max (E) = ^&D{t - to) is 
the source distance giving the maximum flux of CRs of energy E 
at the Earth for a given t and to- Because S for this case follows 
8 oc (t — to) -2 , its value increases or decreases with to(t) for a 
fixed t(to). For to = 10 5 yr the maximum CR source confinement 
time and t = 2 x W 5 yr, the lowest possible age of the unde- 
tected old SNRs adopted in the present study, each point in the plot 
represents the maximum possible anisotropy amplitude at energy 
E due to the old sources in the present model. Calculations using 
other possible values of the source parameters (i.e. to < W 5 yr and 
t > 2 x 10 5 j/r) will give anisotropies that are less than the val- 
ues shown in Fig. 3. In Fig. 3, the solid lines represent the model 
with a — 0.3 and the dotted lines are for a = 0.6; the upper lines 
represent the values calculated at / = 0.3 and the lower lines at 
/ = 0.1. The observed anisotropy data are also shown in order to 
see the extent of the contribution of the old sources. A detailed in- 
spection of the figure shows that the contribution of the undetected 
old SNRs to the CR anisotropy at the Earth cannot be fully ne- 
glected, particularly for the model with a = 0.6. This can be seen 
from the region bounded by the two dotted lines, which contains 
the data points below E < 10 12 eV and above E > 10 14 eV. How- 
ever, here it should be mentioned again that each point given by the 
calculated lines only represents the maximum possible anisotropy 
due to the old sources. Physically, it is highly unrealistic to get such 
a situation at any point of time, and thereby each line represents the 
values observed at the Earth in the whole energy range. However, 
there is always a possibility of observing atleast one of the energy 
points at a given time. Thus, we cannot claim that the undetected 
old sources can fully explain the observed anisotropy data, but we 
can safely conclude that their contribution cannot be completely 
neglected in the study of CR anisotropy at the Earth. 



There also exists a model dependent argument, which supports 
the contribution of the old sources of the observed CRs. For the 
generally excepted particle release time of to = 10 5 yr (which is 
also the value adopted here), all the nearby known SNRs listed in 
Table 1 are quite young (with estimated ages less than 10 yr) ex- 
cept for the Geminga SNR. Thus, they have not yet released the 
CRs into the local ISM. So, it is quite possible that majority of the 
CRs observed at the Earth originate from nearby undetected old 
sources. In fact, the necessity for the presence of such old nearby 
sources has already been disccussed in the spectral study of both 
CR protons (Thoudam 2006) as well as CR electrons (Atoyan et al. 
1995). 
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Figure 3. Maximum CR anisotropy at the Earth due to undetected old SNRs 
with ages t = 2 X 10 yr assuming burst-like injection of particles at to = 
10 5 yr. The calculations are done using Eq. 8 and by taking the source 
distances r = r max (E), where r max (E) = y/ 6D(t — to) is the source 
distance which gives the maximum flux of CRs of energy E at the Earth for 
a given t and to . Each point represents the maximum possible anisotropy <5 
from an old source at a given energy E (see text for details). The solid lines 
are for a = 0.3 and the dotted lines for a = 0.6 in which the upper ones 
represent the values at / = 0.3 and the lower ones at / = 0.1. Data points 
are the same as in Fig. 1. and are shown just to see the extend of the old 
sources contribution. 



5 RESULTS AND DISCUSSIONS 

The analysis presented here involves the study of the effect of dif- 
ferent particle release times from the sources on the CR spectrum 
and the anisotropy below ~ 10 15 eV observed at the Earth. The 
study considers burst-like injection, continuous injection for a fi- 
nite time, injection from a stationary source and energy dependent 
injection as possible forms of particle injection from the sources. 
The study shows that the case of particle injection from a station- 
ary source is the same as that of the continuous injection for a fi- 
nite time period for which ti < t ^ ti (see Case 1 of section 
3.2). Also, for the proton energy range of (0 — 10 15 )eV considered 
here, the energy dependent particle injection model can be repre- 
sented by the burst-like injection of particles at time t = W r 'yr 
except for particles with energies fi^4x 10 14 eV. When applied 
to the nearby known SNRs located within a distance of 1.5kpc, it is 
found that the model with continuous injection of particles gives an 
anisotropy that is too large to explain the observed data at a = 0.6, 
but it predicts a value that is within a factor of 2 for E ^ W 13 eV at 
a — 0.3 and ti ^ 10 4 yr. However, the burst-like injection model 
with a = 0.3 can explain the data more closely at almost two values 
of particle injection time: to ~ 2 x 10 4 j/r, where the anisotropy 
is determined by Monogem for E < 10 13 eV and by the source 
Gl 14.3+0.3 for E > 10 13 eV and ( ~5x 10 4 yr, where Mono- 
gem dominates in the whole energy range below 10 15 eV\ Also, 
the burst-like model with a — 0.6 explains the data effectively 
above 10 12 eV if particle injection takes place at ~ 5 x 10 4 yr, al- 
though below 10 12 eV it predicts an anisotropy that is larger than 
the observed value by a factor of almost 4. For this case, the to- 
tal anisotropy is again determined mainly by the Monogem SNR. 
This result is strongly supported by the close agreement of the ob- 
served phase of the anisotropy with the positions of G 114.3+0.3 
and Monogem in the Galaxy. Therefore, if the observed anisotropy 
is due to the effect of the local sources, the study favours the burst- 
like particle injection model both at a = 0.3 and 0.6. 

The study further shows that the contribution of the undetected 
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old SNRs cannot be fully neglected in the study of CR anisotropy 
at the Earth. Because of their unknown source parameters (such 
as their ages and distances) it is impossible to calculate their ac- 
tual contribution to the observed CRs. However, we can always try 
to make an estimate of their contribution by considering a proper 
model dependent set of source parameters. In the present study, 
using source ages of t = 2 x 10'' yr and particle release time of 
to = 10 5 yr, it has been shown that for the model with a = 0.6 
they can contribute significantly to the local CRs. Calculations us- 
ing other possible values of the source parameters (i.e. to < 10°yr 
and t > 2 x 10 5 yr) will give anisotropies that are less than the val- 
ues obtained here. Thus, Fig. 3 shows only the maximum possible 
anisotropy value at energy E due to the old sources in the present 
model. The study concludes that the contribution of the undetected 
old SNRs to the observed CR anisotropy cannot be simply ignored 
and should be considered in the study of local CRs. 
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